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Abstract. The distribution of times At between coronal mass ejections (CMEs) 
in the Large Angle and Spectrometric Coronagraph (LASCO) CME catalog for the 
years 1996-2001 is examined. The distribution exhibits a power-law tail oc (At) 1 
with an index 7 ~ — 2.36±0.11 for large waiting times (At > lOhours). The power- 
law index of the waiting-time distribution varies with the solar cycle: for the years 
1996-1998 (a period of low activity), the power-law index is 7 ~ — 1.86±0.14, and for 
the years 1999-2001 (a period of higher activity), the index is 7 « -2.98 ±0.20. The 
observed CME waiting-time distribution, and its variation with the cycle, may be 
understood in terms of CMEs occurring as a time-dependent Poisson process. The 
CME waiting-time distribution is compared with that for greater than CI class solar 
flares in the Geostationary Operational Environmental Satellite (GOES) catalog for 
the same years. The flare and CME waiting-time distributions exhibit power-law 
tails with very similar indices and time variation. 



1. Introduction 

Coronal mass ejections (CMEs) are large scale expulsions of plasma 
and magnetic field from the Sun, typically observed in scattered white 
light images of the solar corona and low solar wind. CMEs are closely 
associated with solar flares, explosive events in which magnetic energy 
is released in situ in the solar corona, although the exact relationship 
remains the subject of debate (e.g. Kahler, 1992). 

The statistics of solar flares have been intensively studied. For exam- 
ple, hard X-ray studies suggest that the frequency distribution of flare 
energies (the number of flares per unit energy and per unit time) follows 
a power law, N(E) oc E~ 1 ^ (Crosby, Aschwanden, and Dennis~993). 



The index of the power law does not vary substantially with the solar 
cycle, suggesting that it is of fundamental significance, and it was on 
this basis that the self-organized critical (SOC) or avalanche model for 
flares was proposed (Lu and Hamilton, 1991; for a review of recent 
developments of this model see Charbonneau et al., 2001). 

The statistics of CMEs have also been the subject of scrutiny. Be- 
cause the exact relationship between flares and CMEs remains poorly 
understood, there is fundamental interest in whether flares and CMEs 
exhibit the same statistics. The mass distribution of CMEs appears 
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to follow an exponential distribution ( Jackson and How ard, 19 931), a n d 
this may also be the case for the energy distribution ( ] Jacks'on7T997 ) . If 
this result is confirmed, the difference in energy distributions of flares 
and CMEs may provide a valuable clue to their respective production 
mechanisms. Other statistical properties of CMEs that have been stud- 



ied include the variation of rate with the solar cycle (Webb and Howard, 1994), 
the association with different kinds of soft X-ray flares ( Harrison^ 

the possibility of different CME classes ( Sheeley et al., 1999D, the asso- 

ciation with sigmoidal structures in active re gions (|Ca nfie^^Hucl^om^ami^eyt sov, 2000 ) , 
the association with solar microwave bursts ( Dougherty^rnn^mcl^su^OO^ 
and the correlation of rate of occurrence with properties of originating 
active regions (Falconer, Moore, and Gary, 2002). 

Recently another statistic of flares, the distribution of times At 
between events ('waiting times'), has attracted attention. Boffetta et 
al. (1999) pointed out that the waiting-time distribution (WTD) for 
20 years of Geostationary Operational Environmental (GOES) soft X- 
ray flares exhibits a power-law tail oc (At) 7 for large waiting times 
(greater than a few hours). They argued that the power law indicates 
long correlation times in the flare time series, and attributed particular 
significance to the observed power-law index, 7 ~ —2.4. They also 
argued that the appearance of a power law is inconsistent with the 
SOC models for flares. However, Wheatland and Litvinenko (2002) 
re-examined the soft X-ray flare catalog used by Boffetta et al., and 
found that the power-law index of the tail of the WTD varies with the 
solar cycle, which counts against the power law having a fundamental 
significance. Wheatland and Litvinenko explained the observed WTD 
as arising from a Poisson process with a time- varying rate. Flares occur 
essentially independently, i.e. as a Poisson process in time, but the 
mean rate of the process varies with the solar cycle. Simple theory was 
presented to account for the observed form of the WTD in terms of 
a time-dependent or non-stationary Poisson process. In principle this 
interpretation is consistent with SOC models for flares, and Norman 
et al. (2001) have demonstrated that an avalanche model with time- 
dependent driving can produce a WTD with a power-law tail matching 
the GOES observations. More generally the appearance of power laws 
in the waiting-time statistics of SOC models can be explained in terms 
of correlated driving (Sanchez, Newman and Carrcras, 2002). 

In this paper the waiting-time distribution for coronal mass ejections 
is examined, based on a catalog of events from the Large Angle and 
Spectrometric Coronagraph (LASCO) on the Solar and Heliospheric 
Observer (SOHO) spacecraft. To our knowledge this is the first time 
that waiting times for CMEs have been examined. One motivation 
for this study is to confirm the results presented in Wheatland and 
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Litvinenko (2002) for soft X-ray flare waiting times. More generally, as 
discussed above it is of interest to determine whether flares and CMEs 
share the same statistics, in this case waiting-time statistics. 

The order of presentation of the paper is as follows. In §2.1 the 
catalog on which this study is based is introduced, and the waiting- 
time distribution for all events in the catalog is constructed in § 2.2. 
In § 2.3 the time variation of the WTD is examined, by looking at the 
distribution for the first and second halves of the catalog. In § 3 the 
observed WTDs are compared with those for soft X-ray flares. In § 4 
the observed CME WTDs are explained in terms of simple theory. The 
theory is outlined in §4.1, and is then applied to the CME catalog 
using a piecewise-constant Bayesian decomposition of the CME time 
series (§4.2). Finally, in §5 the results are discussed and conclusions 
are drawn. 



2. LASCO CME waiting-time distributions 



2.1. LASCO CME catalog 

This study is based on the Center for Solar Physics and Space Weather /Naval 
Research Laboratory SOHO/LASCO CME catalog (which is available 
on the web, at http:/ /cdaw. gsfc.nasa.gov/CME_list/). The catalog lists 
all CMEs identified by LASCO operators from the start of SOHO 
observations in January 1996 to the end of 2001, a total of 4645 events. 
For each event the catalog lists the date and time of the first appearance 
of the CME in the LASCO C2 coronagraph field of view, the central 
position angle of the CME, estimates of the CME speed, and a number 
of other details. In this study only the date and time of first appearance 
in the C2 coronagraph is used. 

Figure 1 shows monthly numbers of CMEs in the catalog (solid 
histogram), together with the monthly sunspot numbers (dotted his- 
togram). 1 The CME numbers closely track the sunspot number, as 
pointed out by Web b and Howardl (|1994j) . Figure 1 also indicates by 
vertical dashed lines the start and end of two major data gaps, that 
were due to the temporary loss of the spacecraft in the second half of 
1998, and to a second interruption to the mission late in the year. In 
this study these gaps are taken into account. 



1 The sunspot numbers are available from the National Geophysical Data Center 
at ftp://ftp.ngdc.noaa.gov/STP/SOLAR_DATA. 
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Figure 1. Monthly numbers of CMEs in the LASCO catalog for 1996-2001 (solid 
histogram) and monthly svmspot numbers for the same period (dotted histogram). 
The dashed vertical lines indicate the beginning and end of two major data gaps. 



2.2. WTD FOR 1996-2001 

Figure 2 shows the waiting-time distribution constructed for all events 
in the SOHO/LASCO CME catalog, presented as a log-log plot. The 
distribution exhibits a clear power-law tail oc (At)' 1 for waiting times 
At greater than about 10 hours, and a fit to the distribution (for At > 
10 hours) gives 7 ~ —2.36 ± 0.11. The power-law fit is shown by the 
solid line in the figure. 

Figure 2 confirms that the WTD exhibits similar power-law behavior 
to that observed for the GOES soft X-ray flares (Boffett a et al., 1 999; 



Wheatla nd and Litvinenko, 2002 ) . A detailed comparison with the flare 
data is made in § 3. 

For short waiting times (less than a few hours) the WTD in Figure 2 
exhibits widely varying values in adjacent bins. This effect is most likely 
due to the discrete times between C2 observations. Typical observing 
sequences for the LASCO instrument lead to certain intervals between 
C2 observations being more common, which influences whether CMEs 
that are closely spaced in time are observed. 

2.3. Time variation of WTD 

Figure 1 shows that the observing period (1996-2001) encompasses an 
initial period of solar minimum, followed by the rise and part of the 
maximum of a solar cycle (cycle 23), as defined by sunspot number. To 
examine the secular variation of the WTD, we have divided the dataset 
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Figure 2. Waiting-time distribution for LASCO CMEs for the years 1996-2001. 

into two halves: the period 1996-1998, and the period 1999-2001. The 
first half includes solar minimum and the rise to maximum, and the 
second half is essentially a period of maximum activity. 

Figure 3 shows the WTDs for these two periods. The power-law 
tail of the WTD clearly varies with time: for 1996-1998 the index of 
the power law is 7 f» — 1.86 ± 0.14, and for 1999-2001 the index is 
7 ~ —2.98 ± 0.20. The indices were again determined by fits to the 
observed distributions for At > 10 hours, and the fitting power laws 
are shown in the figure by the solid lines. 

TIME VARIATION OF CME WTD 
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Figure 3. Waiting-time distribution for LASCO CMEs for the years 1996-1998 
(solid), and for the years 1999-2001 (dotted). 
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It is easy to understand, at a qualitative level, the origin of the 
variation shown in Figure 3. During the period of higher activity (1999- 
2001) the mean rate at which CMEs occur is higher, as shown by the 
monthly numbers in Figure 1. The average waiting time is the reciprocal 
of the mean rate, so the average waiting time is shorter during 1999- 
2001, i.e. there are fewer large waiting times. However, this explanation 
does not account for the observed functional forms of the WTDs. A 
detailed quantitative explanation for the observed behavior is outlined 
in §4. 



3. Comparison with flare WTDs 

It is interesting to compare the CME WTDs with those for soft X-ray 
flares for the same periods of time. For the flare data we use the GOES 
soft X-ray flare catalog, which spans the years 1975-2002. 2 The events 
in the catalog were selected from soft X-ray time histories measured by 
the GOES satellites. The soft X-ray background varies by more than an 
order of magnitude over the solar cycle. Consequently small flares are 
not observed at solar maximum owing to the increased background, an 
effect that biases WTDs constructed from all events in the catalog (for a 
more extensive discussion of this point, see Wheatland and Litvinenko, 
2002). To correct for this effect we have restricted the GOES events 
to those exceeding CI class (a peak flux greater than 10 -6 Wm~ 2 
at the spacecraft). This flux is an appropriate threshold because it 
is typical of background values at solar maximum. Note that for the 
CMEs, no similar 'thresholding' is applied, because the detection of 
events is essentially limited by the sensitivity of the instrument, and 
not by selection against a time-varying background. An exception to 
this is that at times of higher activity multiple CMEs occurring nearly 
simultaneously are more common, and in this situation some events 
may be missed due to a 'background' of other events. However, this 
effect only influences the number of short waiting times, and not the 
tail of the WTD. 

Figure 4 shows the WTDs for the CME catalog (solid histogram), 
and for the GOES soft X-ray flares of greater than CI class (dotted 
histogram) for the duration of the CME catalog, i.e. for the years 1996- 
2001. Both distributions exhibit similar power-law tails. Fits to the 
behavior of each distribution for At > 10 hours give power-law indices 
7 « -2.36 ± 0.11 for the CME distribution, and 7 « -2.26 ± 0.11 for 

2 The GOES catalog is available online at 

ftp://ftp.ngdc.noaa.gov/STP/SOLAR_DATA 
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the flare distribution. Hence the two indices agree with one another 
within the stated uncertainties. 



WTDS FOR CMES AND FLARES (1996-2001) 
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Figure 4. Waiting-time distribution for LASCO CMEs for the years 1996-2001 
(solid histogram), and for GOES soft X-ray flares of greater than CI class for the 
same years (dotted histogram). 

As shown in Figure 3, the power-law tail in the CME WTD is 
different for the periods 1996-1998 and 1999-2001. It was shown in 
Wheatland and Litvinenko (2002) that the WTD for the GOES events 
varies with the solar cycle, but how does the time variation of the flare 
distribution compare with that for the CME WTD? Figure 5 shows 
the CME and flare WTDs for the years 1996-1998 (top panel) and for 
the years 1999-2001 (bottom panel). In each panel the CME distribu- 
tion is shown as a solid histogram and the flare distribution is shown 
as a dotted histogram. Power-law fits to the tail of each distribution 
(At > 10 hours) are indicated by solid lines, and by the numerical values 
of indices shown in the figure. For the years 1996-1998 the CME index 
is 7 « -1.86 ± 0.14, and the flare index is 7 « -1.75 ± 0.08. For the 
years 1999-2001 the CME index is 7 « -2.98 ±0.2, and the flare index 
3.04 ± 0.19. Hence we find that the two phenomena exhibit 
power-law indices that vary in essentially the same way with time. The 
observed power-law indices are the same within uncertainties. 
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Figure 5. Top panel: WTDs for LASCO CMEs (solid histogram) and for GOES soft 
X-ray flares (dotted histogram) for the years 1996-1998. Bottom panel: The same, 
for the years 1999-2001. 

4. Explaining the observed WTDs 

4.1. Theory 

The observed waiting-time distributions for LASCO CMEs may be 
understood in terms of CMEs occurring as a time-dependent Poisson 
process, i.e. as a sequence of independent events, with a mean rate that 
varies with time. The necessary theory was presented in Wheatland and 
Litvinenko (2002), and is summarized here. 

A Poisson process is a sequence of events in time such that there 
is a constant probability per unit time A of an event occurring. The 
quantity A is referred to as the rate of the process. The waiting-time 
distribution for a Poisson process is a simple exponential, 

P(At) = Ae~ AA *. (1) 
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A sequence of independent events with a mean rate A = A(t) that 
varies with time t defines a time-dependent Poisson process. Provided 
the rate varies slowly (with respect to the average waiting time), the 
WTD for a time-dependent Poisson process may be written 

P(At) = rjr [ T \(t) 2 e- x( - t)M dt, (2) 
Jo 

where N is the total number of observed events and T is the total 
observing time. Equation © may be recast in the form 

1 r°° 

P(At) = = X 2 f(X)e~ XAt dX, (3) 
A jo 

where /(A) is the time distribution of the rate, i.e. f(X)dX is the fraction 
of time that the rate is in the range (A, A + dX), and where 

- N f°° 

A = - = y o A/(A) dX (4) 

is the average rate. 

For the special case of a piecewise-constant Poisson process, i.e. a 
Poisson process consisting of a series of constant rates Aj for intervals 
t{, Equation (j2j) becomes 

P(At) = ^^e-^. (5) 

i 

To understand the behavior of the tail of the WTD, note that 
Equation © implies that the asymptotic (At — > oo) behavior of P (At) 
depends on the behavior of A 2 /(A) for small values of A. Specifically 
if /(A) has a power-law form /(A) oc A a for A — > 0, then we have 
P(At) oc (At)~( 3+a ) for At -> oo. This result holds for a > -3, and in 
particular includes the case a = 0, when /(A) is a constant for small 
A. Hence a power-law tail in the WTD arises from power-law behavior 
in the time distribution of the rate, for low rates (including the case 
when the distribution is flat for low rates). 

4.2. Comparison of observations with the 
piecewise-constant model 

To apply the theory outlined in §4.1, it is necessary to estimate the 
time variation of the rate of CME occurrence, A = A(t), for events in 
the LASCO catalog. Note that there is no doubt that the CME rate is 
time varying, as evidenced by the secular change in the monthly values 
in Figure 1. 
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To approximate the rate we apply the Bayesian Blocks procedure (Scargle, 1998), 
which takes the time history of CME events and decomposes it into 
a sequence of piecewise-constant rates A, and intervals t{. This pro- 
cedure was previously applied to the GOES flares ( Wheatland, 2 000; 
Wheatland and Litvinenko, 2002 1 . The result of applying the Bayesian 



Blocks procedure to the LASCO catalog is shown in Figure 6. The 
top panel shows the cumulative number of CMEs in the catalog as a 
function of time — the gradient of this graph indicates the mean rate 
of occurrence of CMEs. The two major data gaps (discussed in §2.1) 
are shown by the vertical dashed lines. The bottom panel of Figure 6 
shows the Bayesian rates and intervals. The Bayesian procedure has 
decomposed the six years of observations into 33 intervals during which 
the rate was approximately constant. 
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Figure 6. Top panel: Cumulative number of CMEs in the LASCO catalog versus 
time. Bottom panel: Bayesian estimate of piecewise-constant rates and intervals for 
the LASCO CMEs. In both panels the dashed vertical lines indicate the major data 
gaps. 
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The Bayesian procedure is essentially autonomous, i.e. makes its 
own decisions during the rate decomposition shown in the bottom 
panel of Figure 6. However, a number of free parameters need to be 
specified. First, the method requires observations in discrete timesteps. 
For simplicity we have presented the data in timesteps of five minutes. 
Second, the method involves a prior-odds ratio, which has been set to 
two. This means that the data for a given interval must be considered to 
be twice as likely to be produced by a two-rate process by comparison 
with a single-rate process for the two-rate process to be adopted as the 
model. Third, the minimum number of events in an interval is chosen 
to be three. In practise we find that the rate decomposition is relatively 
insensitive to the choices of these parameters. 

Based on the rates and intervals determined by the Bayesian pro- 
cedure we can construct a model piecewise-constant Poisson WTD, 
using Equation The result is shown in Figure 7 by the solid curve. 
The histogram is the observed WTD for all of the LAS CO events, 
as also shown in Figures 2 and 4. Figure 7 demonstrates that the 
piecewise-constant model provides a reasonable representation of the 
observed WTD, including reproducing the extended tail of the distri- 
bution. It should be noted that the piecewise-constant model involves 
a fairly crude approximation — the true rate of occurrence of CMEs 
is expected to be continuously varying. Despite this limitation, and 
the approximate method of rate determination, the piecewise-constant 
Poisson model is seen to reproduce the observed form of the WTD. 



WTD FOR LASCO CMES (1996-2001) AND MODEL 
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Figure 7. WTD for the LASCO CMEs for 1996-2001 (histogram), together with 
the WTD for the piecewise-constant Poisson model, with rates and intervals taken 
from the Bayesian estimates in the bottom panel of Figure 6. 
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It is straightforward to construct the piecewise-constant model (jSJ) 
for the periods 1996-1998 and 1999-2001 separately, based on the Bayesian 
estimates in the lower panel of Figure 6. The result is shown in Figure 8. 
The histograms show the WTDs for the two periods, as also shown in 
Figure 3. The curves show the piecewise-constant models for each pe- 
riod. The piecewise-constant models are seen to qualitatively reproduce 
the observed WTDs. Once again, it should be noted that although the 
models are crude, the observations are reproduced reasonably well. 




Figure 8. Observed LASCO CME WTDs for 1996-1998 and for 1999-2001 (solid 
histograms), together with piecewise-constant Poisson models (solid curves). 



5. Conclusions 



In this paper the waiting-time distribution for coronal mass ejections is 
examined, based on the LASCO CME catalog, which spans the years 
1996-2001. The main results are first that the WTD for the entire 
catalog exhibits a power-law tail (for waiting times At larger than 
about 10 hours), i.e. the distribution has the form cx (At) 7 , with a 
power-law index 7 ~ —2.36 ± 0.11. Second, the power-law index of the 
tail varies with time, so that the power law is significantly steeper at 
times of higher activity, i.e. when CMEs are more frequent. Third, the 
power laws in CME WTDs are the same as those found in the WTDs of 
GOES soft X-ray flares, for the same years (i.e. flares and CMEs exhibit 
the same power-law tails in their WTDs, and the same time variation 
of those tails). Finally, the observed WTDs can be explained in terms 



ms.tex; 2/02/2008; 3:13; p. 12 



CME WAITING TIMES 



13 



of CMEs occurring as a Poisson process in time, with a time-varying 
rate. 

The theory presented in §4.1 describes the observed WTDs as aris- 
ing from a Poisson process with a time-dependent rate A = X(t). This 
means that the theory assumes that CMEs occur independently of one 
another, and with a probability per unit time that varies with time. It 
is necessary to explain both of these aspects of the model. Concerning 
the independence (randomness) of CME occurrence, the first point to 
note is that consecutive CMEs in the catalog may often originate from 
different regions on the Sun, which are expected to produce CMEs 
independently. However, at times of solar minimum the majority of 
CMEs may originate from one or a few active regions at any given 
time. Although the detailed theory of the CME mechanism remains 
to be worked out, numerical models often involve complex magnetic 
configurations that are presumably reached via complicated processes 
of coronal magnetic field evolution (e.g. Amari et al., 2000). Such pro- 
cesses may explain the origin of the apparent independence of successive 
CMEs from a common region on the Sun. Regarding the time-varying 
rate of CME occurrence, on short timescales the CME rate is expected 
to vary due to the appearance and disappearance of active regions with 
suitably complex magnetic configurations. On longer timescales the 
rate is expected to vary with the solar cycle due to the cyclic variation 
in the number of active regions (see Figure 1). 

As noted in the introduction, the poor understanding of the re- 
lationship between flares and CMEs means that there is interest in 
whether flares and CMEs exhibit the same statistics. In § 3 of this 
paper it is shown that flare and CME waiting-time statistics are very 
similar, and in particular the two phenomena exhibit the same power- 
law tails in their WTDs. The theory in §4.1 of this paper shows that, 
assuming the occurrence of flares and CMEs can be understood as a 
time-dependent Poisson process, the observed WTD depends only on 
the time distribution of the rate of events. In this case the appearance 
of similar WTDs implies that both flares and CMEs occur with a rate 
that varies with time in a similar way. At face value this result suggests 
a common origin for flares and CMEs. However, it should be noted that 
there are pieces of observational evidence suggesting distinct origins for 
flares and CMEs. For example, at solar maximum a significant number 
of CMEs originate from very high solar latitudes, well away from active 
regions ( Hundh ausen, 19 99 ) . 

The appearance of different power-law tails in the CME WTDs was 
explained in §4.1 in terms of power-law like behavior in the time dis- 
tribution of occurrence rate at low rates. At times of high activity the 
observed WTD has an index close to 7 = —3, which implies that the at 
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these times the time-distribution of rates is approximately flat at low 
rates. For times of lower activity the WTD has a shallower power law, 
implying that at these times there is a decreasing power law in the time 
distribution of rates at low rates. This behavior is consistent with that 
previously found for flares (Wheatland and Litvinenko, 2002). The ori- 
gin of this interesting difference in CME (flare) behavior between dif- 
ferent phases of the cycle remains an open question. 

The absence of a fixed power law in the tail of the CME and flare 
WTDs suggests that the power law does not have fundamental signif- 
icance, such as that attributed to the power law in the flare energy 
distribution. The theory outlined in §4.1 explains the observed power 
laws, and suggests that they will appear generically for phenomena 
involving independent events with a time-varying rate. However, the 
power laws do have some significance, in that they contain important 
information about the rates at which CMEs and flares occur. In partic- 
ular, the implication of power laws in the underlying rate distributions 
may be a valuable clue for understanding the mechanisms of energy 
storage and release in flares. 
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